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Abstract 

We analyze long range rapidity correlations observed in the STAR experiment at RHIC. Our goal 
is to extract properties of the two particle correlation matrix, accounting for the analysis method 
of the STAR experiment. We find a surprisingly large correlation strength for central collisions 
of gold nuclei at highest RHIC energies. We argue that such correlations cannot be the result of 
impact parameter fluctuations. 
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1. Introduction 

The STAR experiment at RHIC has reported preliminary observations of forward backward 
correlations as a function of the centrality of the collision The reported forward-backward 

correlation shows a rapid increase in strength as a function of centrality, and appears to have a 
strength which cannot be explained by a superposition of pp interactions. It has been suggested 
that such correlations are a consequence of the Color Glass Condensate induced Glasma produced 



in the early stages of heavy ion collisions 0, i, i, S, 0, S, 0, EH, E, El, Q, [H, EE! 



^ ~ ~ fz} — — ^ i^-^j t^^j I ~ ij I ij I 17 I I I I ' i-hU-i- 

An alternative but quite similar explanation is provided by parton percolation [2l|, 
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long range correlations may be the origin of the ridge phenomenon measured at RHIC, ^ 

23,13, 0,113, |3l| that may also be explained as arising from color electric and magnetic flux tubes 
originating in the Glasma [32l . 

The results for the forward-backward correlations strength measured in the STAR experiment 
are reported as a function of centrality bin. One might be worried that if one averages over all 
the events in a centrality bin, correlations would be generated by the different impact parameters 
(or numbers of wounded nucleons) possible within such a bin; this is essentially the mechanism 



proposed in 



37 



We expect the charged multiplicity to be strongly correlated with 



impact parameter and if the impact parameter itself can have significant variation within a fixed 
centrality bin, then spurious correlations whose only origin is the geometry of the collision would 
be generated. 

The STAR analysis is however more sophisticated. For each event one measures also a reference 
multiplicity in a relatively wide central rapidity interval. Correlations and fluctuations of the 
forward and backward multiplicities are then measured separately for each reference multiplicity. 
The correlations and fluctuations reported are average values of these measurements over each 
centrality bin. One might expect that such a procedure is relatively insensitive to impact parameter 
correlations. 
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One purpose of this paper is to show how to extract values of the multipUcity correlation 
function including the constraints actually applied by the the STAR experiment. After setting 
up some notations in sec. [2] we argue in sec. [3] that, within a Gaussian approximation for the 
fluctuations, there is a bound on the ratio of forward backward correlated fluctuations to forward 
fluctuations, 

^ _ {NfNb)n^ - {Nf)n^{Nb)n^ 

Here, is the multiplicity measured in some forward bin of rapidity, and A^b the value in a 
backwards bin, and the subscript Ar indicates that the averages are taken for fixed values of 
the reference multiplicity. In the STAR experiment the pseudorapidity windows are at symmetric 
values of pseudorapidity around r/ = 0. When the forward and backward rapidity windows are 
smaller than the reference window and farther from each other than from the reference window, 
we argue on quite general grounds that b < 1/2. This is perhaps violated by non-linear effects 
that are of order 1/(Af), effects that could be several percent for the most central collisions. If 
b = 1/2, the multiplicities are maximally correlated. This bound is derived making, besides this 
Gaussian approximation, only minimal assumptions about the correlations and does not depend 
on the detailed mechanism for generating them. The STAR data is consistent with this bound 
except for the most central events of Au-Au collisions at the highest beam energy, where there are 
small violations. We then analyse the effect of the impact parameter induced correlations of the 
type suggested in Ref. 36, 37, H, 39] in sec. [Hand find that they are not sufficient to explain the 



large correlation seen in the STAR data. In sec. [5] we then propose a simple parametrization to 
add intrinsic long range correlations to the impact parameter fluctuations. 

Our results are: The correlations measured in STAR do not appear to be entirely associated 
with impact parameter fluctuations. The measured value of b is near its maximally allowed value 
for central Au-Au collisions at RHIC energy. Indeed for the highest energy and most central Au-Au 
collisions it is larger than the bound b < 1/2. We find that a crucial requirement to obtain a b 
close to this limit are large fluctuations in the reference multiplicity. 

2. The rapidity correlation function 

The quantity that we wish to extract information from experiment is the rapidity correlation 
function. We define it as 

/dNdN\ /diV\/dAr\ 

This can be decomposed into a local piece that corresponds to Poissonian fluctuations in the particle 
number and a long range correlaion 

C(,-,') = %-,')(f) + A-h-,/)(^^)(f;j). (3) 



The 6(7] — 7]') piece is often (e.g. 40|, l4]|) absorbed into the definition of the correlation function 
itself (the l.h.s. of eq. ([2])), but it will be convenient for the following discussion to keep it on the 
r.h.s. as part of C(r/ — r/'). 

Some information about the function K can be found by measuring local multiplicity fluctua- 
tions. Integrating eq. ([3]) twice over a pseudorapidity interval Arj centered at some rapidity tj and 
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assuming that Ar] is sufficiently small that the pseudorapidity dependence of the multiplicity in 
this interval can be neglected, we obtain 



(n2) - {nf = (n) + 




(4) 
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This result is identical to the result for obtained a negative binomial distribution, and we can 
identify the A;-parameter as 



(The negative binomial distribution has further consequences for higher moments of the multiplicity 
distribution as well.) This relation can be used to gain information on the rapidity correlation 
for small rapidity differences from measurements of the multiplicity distribution in different Atj- 



Of central interest in this paper is the behavior of the correlation function C at large rapidifr 
differences. In particular, we shall concentrate on the case studied by the STAR collaboration 



This analysis measures the correlation between charged particle multiplicities in forward ("F") and 
backward ("B") rapidities. Specifically, pseudo-rapidity windows of width 5 = 0.2, are situated 
symmetrically around midrapidity within the STAR TPC acceptance \r]\ < 1. The data is di- 
vided into 10% centrality bins and the correlation coefficient between the forward and backward 
multiplicities is measured within the centrality bin. In order to eliminate the effect of centrality 
fluctuations within one bin, there is an additional crucial twist in the analysis. One also measures 
the multiplicity in a third reference ("R") rapidity interval of width (5r = 1.0 that does not overlap 
with the "F" and "B" windows. The variances and the covariance of the "F" and "B" windows are 
then determined separately for each reference multiplicity N^, the motivation being that a fixed 
selects, with a good accuracy, events with a fixed impact parameter. Although the data is 
presented for rapidity intervals in different locations, it is good to keep in mind a typical config- 
uration, with rapidity windows -1 < r/ < -0.8 ("B"), -0.5 < r] < 0.5 ("R") and 0.8 < r/ < 1 
("F"). 

3. Gaussian approximation 

As we have discussed, in order to analyze the STAR results, we shall have to study the case 
of three correlated multiplicities, N-p, Nb and Nji. We shall first compute such fiuctuations in 
a Gaussian approximation for the probability distribution of these three variables within all the 
events in a centrality bin. This should be a good approximation when considering fluctuations 
around an average multiplicity, so long as that average multiplicity is large. This is a consequence 
of the central limit theorem, and in sec. [H we shall see explicitly how it follows directly from a 
stationary phase approximation to expressions for the multiplicity fluctuations in a simple model. 
The corrections to this approximations should be of order 1/A^, where is a multiplicity in a 
bin. For the STAR experiment, and central events, this should be valid to a few percent, since 
the typical value of Np^B ~ 100 for the forward or backward multiplicity. However, in events with 
centrality around 50%, Np^s ~ 20, so one can expect 10-20% corrections. 



We assume that there is some requirement on the centrality of the collision that flxes the average 
values (iVp), (Nb), and (Ar). We define Au = Nu - (Nu), where C/ is F, B or R, so that the 




(5) 
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normalized probability distribution is 



PiNF,NB;N^) 



(27r)3/2 det E 



exp 



(6) 



In this equation 
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{AuAv) = / dr? / dr]' C{r] - ??') 



U,V = F, B, R. 



(7) 



In all that follows, we will assume that the forward and backward windows are chosen to be 
centered at rapidity values symmetric around r] = 0, and that the reference multiplicity is chosen 
by summing over rapidity values that are also symmetrically displaced. Since ^[/y = T,vu by 
construction this means that 



(8) 
(9) 



Other than the Gaussian form and the symmetry above, we shall in this section make no assump- 
tions about the physical origin or strength of the correlation; except that it is a decreasing function 
of the rapidity separation between the measured multiplicities. We also introduce the correlation 
coefficient between two multiplicities with the conventional definition 



Ruv 



uu^vv 



(10) 



The correlation coefficient is mathematically restricted (due to the Cauchy-Schwarz inequality) to 
values between —1 and 1. Without fixing the reference multiplicity, the covariance Df^f, variance 
Dff and correlation coefficient b used in [H, [3] would be the same as 0"^^ a"^^ and R^^- We 
shall however reserve the notation Df^f Djf and b to the quantities at fixed that are actually 
measured. The probability distribution ([6]) involves S^y, the inverse of the correlation matrix S. 
As a 3 X 3 matrix S is easily inverted to give 
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When the reference multiplicity is not measured, the probability distribution can be reduced 
to the double distribution for Ap, Ab by integrating over Ar, 



P(Af,Ab) = y"dARP(AF,AB,AR) 



(27r)3/2 det Sa 



exp 



U,V = F,B 



(12) 



where S2 is the inverse of the 2x2 correlation matrix 



^FF ^FB 



(13) 



In our case, however, we want the probabihty distribution of Af, Ab for a fixed Ar, not integrated 
over all values of Ar. In this case the (conditional on Ar) probability distribution for Ap, Ab is 
given by 



P(Af,Ab|Aj 



P(Af,Ab,Ar) 



P(Aj 



(14) 
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27r det S 
U,V 



exp 
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(Srr)- 



"^R 



F,B. 



The presence of mixed terms in AfAr,AbAr means that the expectation values of Af,Ab are 
shifted from zero and proportional to Ar. This is natural, since if one has an event with exception- 
ally large A'^r one also expects a large A^f and A^b because the multiplicities are correlated. The 
term in A^ assures the correct normalization J dAFdAB-PR(AF, Ab) = 1. The fluctuations and 
correlations of Af, Ab are determined by the coefficients of the quadratic part. They are described 
by the reduced correlation matrix S, obtained by taking the F,B-elements of from eq. (Ilip 
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To get the correlation matrix of Af , Ab for a fixed Ar we must then again invert this matrix to 
get 
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From this we can read off correlation coefficients that correspond to the measured observables 
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The shifts in the expectation values of the A's for a fixed reference multiplicity can be computed 
in terms of the matrix elements of T, as 



(Af, 



B)Nn 



-A 



R 



(20) 



The information in the correlation between the reference multiplicity and the forward and backward 
multiplicity, as well as the fluctuations in the reference multiplicity itself contain some non-trivial 
information about the correlations. 
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Equation (|19|) is the central result of this section. It describes how measurement of the forward- 
backward correlation is modified by the fact that the multiplicities are observed for a fixed reference 
multiplicity, which in turn is correlated with the forward and backward multiplicities. We shall 
now turn to analysing its consequences using some very general assumptions on the two particle 
correlation ([2]). The corrections caused by the fixed reference multiplicity have a very intuitive 
meaning. For a fixed reference multiplicity the Ap^B fluctuate less, SpF < fpp- This is due 
to the correlation between the reference and F, B multiplicities and the effect goes away when 
-^FR ~^ 0- The modification also decreases when the reference multiplicity has larger (uncorrelated) 
fluctuations (T^„ . 

In the limit of very strong correlations, the correlation coefficients R are all very large. This is 
the case in heavy ion collisions when the centrality bins are taken to be very wide so that fluctuations 
and correlations are dominated by impact parameter fluctuations within the centrality bin. Let us 
first assume, for simplicity, that all the correlation coefficients are equal, R^p = R^^ = R^^ = R. 
Now the conditional correlation coefficient (|19p is given by 

R — R^ R , , 

Remembering that R <1 this lead to the upper limit b < 1/2. 

Consider now the typical experimental situation mentioned above, with rapidity windows — 1 < 
r] < -0.8 ("B"), -0.5 < ?7 < 0.5 ("R") and 0.8 < r/ < 1 ("F"). We shah only assume that the 
correlation function C{ri,r]') is boost-invariant and symmetric, i.e. a function of \r] — rj'\ only. Any 
reasonable correlation function C{\r] — r]'\) will decrease as a function of \ri — r]'\. These very general 
assumptions inserted into eq. d?]) lead to simple inequalities for the correlation matrix: 

4"L, < ^-FF = ^^BB '^R > (22) 

Because in our typical configuration the F, B- windows are closer to the reference rapidity window 
than to each other, the average value of the correlation function between the F, B windows and 
the reference window is also larger. This leads to the estimate 

^-L>^-L- (23) 

Together eqs. ((22]) and (f23]) imply that R^^ = R^^ > R^p. Because the constrained correlation 
coefficient b in eq. (fT^ is a monotonously decreasing function of R^^ = R^^^, we then get the upper 
limit ^ ^ 

^ ^ ^BF - -^FR < ^BF - -^BF ^ ^BF < (3 5 > ^_ ^24) 

1 ~ -^FR 1 ~ -^BF 1 + -^BF 

The preliminary STAR results [l|, 0,143] violate this bound derived from very general assumptions, 
and at the moment we see no compelling explanation for this discrepancy. Note also that for the 
configurations where the F, B- windows are close to each other, their distance from each other can 
be less than that to the reference rapidity region. When F, B windows are closer than 0.6 units 
in rapidity, the reference window is determined as 0.5 < \rj\ < 1.0 and eq. ()23p no longer strictly 
applies. In this region the observed value of b does indeed tend to be slightly larger than for large 
rapidity separations. 
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4. Poissonian and impact parameter fluctuation model 



.1. The model 



For an illustration let us consider the toy model introduced in Ref. 36|], where the authors 
claim that the impact parameter fluctuations within the centrality bin would explain most of the 
STAR correlation measurement. The model includes a short range Poisson correlation and the 
fluctuation of the collision geometry via a Monte Carlo Glauber (MCG) model. In practice the 
charged particles produced in a nucleus-nucleus collision are assumed to be a superposition of 
particle production from each participant nucleon: 



iVch = X]ni, (25) 



i=l 

where A'p is the number of participant nucleons. The variables rij are independent and distributed 
according to a Poisson distribution with mean n. These particles are then distributed in pseudo- 
rapidity with a gaussian distribution so that 

4-2. Gaussian approximation 

Knowing that the only correlations in this toy model are the local one leading to the Poisson 
distribution and an infinite range one from the fluctuating number of participant nucleons, one can 
immediately write dowr0 the correlation function of the model: 

c(.,'.^.,-..(f).^(f>(f), 

where dp is the variance of the number of participants in the given centrality bin. It is easy to see 
that this correlation function, when integrated over some intervals in the pseudorapidities i] and 
rj' , reproduces the variance of A^ch defined by eq. (p5]) . 
Integrating this over the rapidity windows gives then 

2 



^PB = 777T2(^f)(^b) (29) 



<j- 

and similarly for the other interval combinations. This gives correlation coefficients 

^BF = ^ , (30) 

^ + (iVp> y ^ + (7Vb> 



^ We are neglecting quantities that are higher order in the multiphcity fluctuations. For example one is approxi- 
mating (n^) ~ {ni)^ in the second term, since it is already proportional to the relative variance a^/ (Np) 
be assumed to be small. 



and correspondigly for R^^ and R^^ . 

To get a rough estimate of the numbers we shall neglect the r/-dependence in eq. ()26p and 
approximate n = n/ y^Tno^ ~ 2 (this is the charged multiplicity per unit rapidity and per partici- 
pant). Taking the rapidity intervals 6 and 5r for the different windows we then get 

^BF = ^ 1 (31) 
-^FR = ; , / =^ (32) 

1 _| 3_ /i _i 1 

ujpiiS y ' uipfiSji 



where we have, following 361], denoted cOp = a^/ {Np). We then take from [Sa] the values Up « 2.5 
for the 0-10% bin and Up w 1.3 for a midcentral 40-50%. bin. This gives R^^ 0.5 and R^^ 0.65 
for the more central bin and Rr,„ « 0.34 and R„„ 0.5. These values of are the correlation 



coefficients quoted in [361]. As we have discussed, however, they do not correspond to the actual 
measured quantity. 

Only at this point need we evoke the Gaussian approximation for the probability distribution. 
If we assume that the probability distribution is Gaussian, we can use eq. (fT9l) to evaluate the 
b parameter from the correlation coefficients ()3ip and (j32p of the model. Putting the numbers 
estimated above into eq. ()19p gives 6 ~ 0.14 for the central and 6 ~ 0.12 for the midcentral bin. 



These values are clearly far from the experimental result. 

4-3. Beyond the Gaussian approximation 

Although it seems unlikely, based on the result we just obtained, that this simple toy model 
would come near to explaining the experimental data, we shall still continue analysing it further 
in order to understand the effect of the Gaussian approximation. For this it is instructive to start 
over from the definition of the model as a probability distribution. We shall again neglect the 
ry-dependence of the single particle spectrum to avoid encumbering our notations. The model is 
defined by a fluctuating number of participant nucleons A'p, which gives the parameter for the 
Poisson-distributions of the measured multiplicities Np, Nb, N^. We can therefore write down 
the probability distribution that characterizes the model as follows 



P{Np, Nf, Nb,Nb) = PucciNp) 



m J \ m ) \ ^vr! 

(33) ' 

Here -Pmcg(A'p) is the probability distribution of events with different numbers of participant 
nucleons, presumably to be obtained from a Monte Carlo Glauber calculation. The first step in 
our approximation is to replace the Poissonian distributions by Gaussians in A'^p, N^, A^r, which 
should be pretty safe as long as these multiplicities are all large enougbH. We then obtain 



{Nf - 5nNpf _ [Nb - 5nNpf _ {N^ - SnnNpf 
25nNp 26hNp 2(^RnA'p 

(34)- 



■^We are simultaneously approximating the discrete variables A^f, Nb, N-r, by continuous ones 
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Figure 1: The correlation coefficient b resulting from a direct Monte Carlo-Glauber evaluation of the probability 
distribution (|33p . Also shown is the result of the approximation where the Poisson distributions are replaced by 
Gaussians, but the distribution of A'^p values is the same MCG one, i.e. eq. (I34p . 



The next step in our approximation introduces a larger error. In the experimental analysis the 
centrality class is defined by events where the reference multiplicity A'r falls between a given lower 
and upper bound. Since the number of participants is not a directly experimentally observable 
quantity, it cannot be used as a selection criterion. We shall, however, assume that the centrality 
bin is defined by a Gaussian distribution of A^p with some variance dp: 



PucciNp) 



1 



exp 



(35) 



The remaining approximation is to replace A'^p by (iVp) in the variances (but not the means!) of 
the Gaussians in Ny, Nb, A^rH- This has the result of turning our probability distribution into a 
Gaussian in all four of its variables: 



PiNp,NF,NB,NR) 



exp 



(Ap-(Ap))2 



{Np-ShNpf _ {NB-ShNpf 



2Sh{Np) 



25n(Np) 



25Rn(Ap) 



(2^)2(rx(Arp))3/25^ap 



(36) 



The remaining step is now to integrate eq. (j36p over A^p. Note that this integration is now done 
for fixed Np, Nb, iVpu- If, as is typically the case for a realistic centrality bin, cTp is relatively large 



This is equivalent to the approximation (n^) ~ (rii) in the variance of dN/drj mentioned in an earlier footnote. 
''Note that since at no point do the experimentalists measure directly the impact parameter or Ap, the correct 
procedure to get a distribution that corresponds to the measurement is to ffist integrate P{Np, Ny, Nb, A^r) over Ap 
to get the probability distribution P{Nf , Nb, N-r). From this one can then restrict to the conditional distribution 
P(Af, A^BlA^Rfixed) = P{Nf, Nb, N-i{,) / P{Nn) and compute Dff and -D/b. This is not the same thing as computing 
Dff and Dft for a fixed impact parameter or Ap and then averaging the results over A^p; which would imply that 
Ap is actually measured event-by event. 
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and 5r ^ 6, this integral is dominated by the reference multipHcity term {N^ — 5jihNp)^. This 
means that selecting events with a fixed selects events with A'^p close to N^/{6jih), not the 
typical value (A'^p) of the original set of events. This causes the probability distribution of A'^p and 
A^B to be peaked around {6/5ji)Nii and thus induces a correlation between Np, N-q, N^. An easy 
way to perform the A^p-integration is to realize that the result is a Gaussian in Np, N-q, Nji. To 
obtain the coefficients of the Gaussian in A'^p, A'b, A^r it is then sufficient to evaluate the expectation 
values and the covariance matrix, which can also be done by integrating first over the multiplicities 
and only then over A'^p. In any case, from eq. ()36p one obtains 

{Nf) = {NB)=6n{Np) (37) 
(Nn) = 6Kn{Np) (38) 

2 

(Arp2) - (ArF)2 = 5fi{Np) + {5hfal = {Np) + j^{NFf (39) 



2 



(ATpATb) - (Arp)(ArB) = {5fif al = j^{Nf){Nb) (40) 

and correspondingly for the other combinations of F, B, R. We have how rederived the formulas 
used in eqs. (|28|) and (|29|) to evaluate the result of the model of in the Gaussian approximation 
previously. 

In order to evaluate the importance of the different approximations we have also used a simple 
Monte Carlo Glauber implementation to evaluate the probability distribution (I33p . The details 
of our simple Monte Carlo Glauber model are described in Appendix |X] and the result of the 
computation is shown in fig. [H It turns out that the result is very close to b = d/{6 + 5r,) = 1/6. 
This corresponds to the limit cjp ^ cxd in the results eqs. (j3ip and (j32p of the toy model. This 
is easy to understand a posteriori. In the Gaussian approximation we assumed that one is only 
looking at events in a centrality bin that is defined by a relatively compact distribution in A'p. 
This is in fact not the case; the centrality bin is defined using the reference multiplicity, and the 
distribution of impact parameters -Pmcg(-^p) in eq. (f33]l includes the whole "min bias" distribution 
of impact parameters, and is therefore very wide. The fact that one must take the limit dp — > cxd 
means, in the language of Sec. El that the correlation coefficients are all approaching unity. They 
do so, however, in such a way that b approaches (5/(5 + (5r). This is in some sense the natural upper 
limit for 6, if the correlation function eq. ([3]) consists of only a Poissonian and long range piece; 
and the most realistic way of modifying it is to introduce a short range correlation which increases 
the fluctuations in the reference multiplicity. We shall now turn to a more general parametrization 
where this can be seen more explicitly. 



5. The effect of intrinsic correlations 

Let us then construct a simple parametrization of long range correlation effects in addition 
to the impact parameter fluctuations. As we saw in the previous section, assuming a Gaussian 
distribution for the impact parameter fluctuations (i.e. neglecting the fact that the centrality 
selection is done using Nji) is not a very good approximation. We would therefore like to construct 
a model where the impact parameter fluctuations are parametrized by an A'^p drawn from a Monte 
Carlo Glauber calculation, and particle production for a fixed A'^p then includes physical long range 
rapidity correlations. We emphasize that we are not assuming that the physical mechanism of 
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Figure 2: The constrained correlation coefficient & as a function of A^r for K — 0.001 and different values of a. 



particle production could be decomposed into independent production from participant nucleons. 
We use the quantity A^p as a convenient proxy for the impact parameter dependent overlap area of 
the colliding nuclei and its fluctuations. The main reason for using A^p instead of the overlap area 
is its easy implementation in our simple Monte Carlo Glauber code. 

We will assume that the probability distribution for a fixed A^p is Gaussian: 



1 



(27r)3/2 detS(iVp) 



exp 



(41) 



where now Au = Nu — Ni/{Np) (we denote the expectation value for a fixed A'p by Nu{Np) to 
separate it from the measured expectation value (Nu) which is averaged over some range of impact 
parameters). The expectation values for a fixed A'^p are 



Nf = Nb = N = 6hNp 
Nk = 6KnNp 



(42) 
(43) 



and also the correlation matrix for fixed impact parameters depends on A'p through the expectation 
values Nu as 



/ Nf + {K + a)N^ 



S(iVp) 



1 



detS 



V 



KN 
KNNr 



N + {K + a)N 
KNNk 



KNNyi 



\ 



(44) 



ATr + (K + a)NB^ I 



Here we have introduced two additional parameters: a describes the increased fluctuations of 
the multiplicity due to a short range rapidity correlation, and K represents a dynamical long 
range fluctuation that increases both the local fluctuations SpF, • • • and generates a correlation 
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Figure 3: The constrained variance Dff and covariance Dbf with the same values for a and K — 0.001 as in fig. E) 



between the different rapidity windows SpR, 
correspond to a parametrization 



In terms of the rapidity correlation function these 



= 5{ri-r]') 



/dN\ 



K + a9(Y-\r]-r]'\) 



/dN\ /dN\ , , 



where y is a characteristic scale of short range rapidity correlations that we assume is larger than 
the size of the rapidity windows but less than the separation between different rapidity window^. 
Our parametrization reduces to the model of Ref. [36t | in the limit a = if = 0. The natural scales 
at which the values of a and K can vary can be estimated from the relation to the parameter k of 
the negative binomial distribution of multiplicities. It can be seen from eq. ([5]) that for relatively 
small rapidity intervals Xjk = K (recall that this is true for a fixed impact parameter, or fixed 
A'p in our parametrization). For central collisions the PHENIX experiment [4^] k ~ 690 for central 
collisions, so typically we would expect a, if ~ 0.001. 

We can then evaluate the correlation using a simple Monte Carlo Glauber implementation 
described in Appendix [XI Our results are summarized in figs. El [3l [H Figure [2] shows the 
dependence of h on a; the strength of the short range rapidity correlation. We see that increasing 
a increases the correlation h. This is to be understood in the following way: the main effect of 
increased short range correlations a is to increase the fluctuations (uncorrelated with A'^p, A^b) in 
Ar. These increased fluctuations then decrease the correlation between the reference and the F,B 
multiplicities. As can be seen from eq. (|19|) . this has the effect of increasing h. For a very large 
a. ~ 0.05 one can even reach h ~ 0.5, but as can be seen from fig. [3] the corresponding values 
of Dff,Dfi, become much larger than the experimental values. Increasing a increases the ratio 
b = Dfb/Dff closer to the experimental result, but it also increases the fluctuations in Np and Ab 
separately beyond what is observed experimentally. 

Figured] shows the dependence of the correlation on the parameter K describing the dynamical 
long distance rapidity correlations. It is seen that both Djf and Dfh and therefore also b are. 



^ We are neglecting the fact that in the actual STAR results the reference window is often wider than its distance 
from the F,B windows. Taking this into account folly would unneccessarily copmplicate the parametrization (|44|l 



without changing our results much. 
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Figure 5: The Monte Carlo result for b with a and K depending on centrality as K — 1/Np and a ~ with a 

coefficient shown in the figure. 
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at the statistical accuracy of our Monte Carlo calculation, independent of K. The most natural 
explanation for this is that, as observed previously, the impact parameter fluctuations already 
generate a correlation that corresponds to the limit dp ^ cxd, or a maximal FB correlation. Adding 
a small physical dynamical correlation described by a reasonable K is negligible compared to this. 
This leads to a significant observation concerning the experimental analysis technique, namely that 
it would be better if the centrality selection could be done, into as narrow bins as possible, using 
another observable as independent as possible from N-^. An example of this effect is the STAR 
observation (fig. 1(a) of Ref. [H) that determining the centrality using the ZDC instead of A^r 
leads to a smaller measured h. 

In our simple parametrization we have so far completely neglected the dependence of the pa- 
rameters on centrality. The charged multiplicity per participant n varies among centralities and, 
more importantly, the dynamical correlation strengths a, K should to a first approximation scale 
with the inverse transverse overlap area, or ~ 1/A^p. Figure [5] shows the result for b when this kind 
of a scaling is taken into account. 

The purpose of this paper is to stress the qualitative effect of taking into account the correlation 
with the reference multiplicity, not to perform a detailed fit to experimental data. A further 
finetuning of the centrality dependence of the parameters would not change the features we are 
addressing here. 

6. Conclusions 

The STAR measurements of long range rapidity correlations point to a very intriguing picture 
of strong correlations from the initial strong color fields in the initial stages of the collision. In 
spite of the apparent simplicity of the experimental observable (counting charged particles in a 
relatively large region of the detector), the measurements turn out to be challenging to interpret. 
The experimental analysis is done by treating separately events with different reference multiplic- 
ities. Thus for a consistent treatment one must consider on equal footing also the correlations 
with the reference rapidity window. This turns the problem from a 2-variable into a 3-variable 
correlation, which has not always been fully appreciated in the literature. The values quoted by 
the STAR collaboration are for a forward-backward correlation for a fixed reference multiplicity. 
Since the result shows that for central collisions there is a strong correlation between the forward 
and backward multiplicities it would be unphysical to neglect the correlation with the reference 
multiplicity. 

We have in this paper discussed long range rapidity correlations in the charged particle mul- 
tiplicity in terms of only very general assumptions on multiplicity correlations. In a Gaussian 
approximation for the probability distribution of events in a centrality class, we find an upper 
limit 6 < 1/2 for the conditional correlation when (5 < 5r and the F and B windows are far 
from each other in rapidity. We then give up the Gaussian approximation and construct a simple 
parametrization of the long and short range correlation, including the effects of impact parameter 
fluctuations with a simple Monte Carlo Glauber model. We show that impact parameter fluctu- 
ations alone are not sufficient to explain the observed data. Because of the correlation with the 
reference multiplicity, the measured conditional correlation coefficient h turns out not to be very 
sensitive to the strength of the dynamical long range correlation. It does depend strongly on the 
short range rapidity correlation through its effect on the fluctuations of the reference multiplicity. 
Nevertheless we do not find a parametrization that would agree with the large h reported by the 
STAR experiment. More experimental data would be welcome to disentangle the interplay between 
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impact parameter, short range correlations leading to increased fluctuations and genuine long range 
correlations. Having a wide enough rapidity between the measured windows to actually see the 
decrease of the long range correlation (as opposed to the geometrical ones that are truly infinite 
range), would be useful but might not be possible within the STAR TPC. The dependence of the 
correlation on the sizes of the rapidity windows {5 and 5^) could be revealing. Preliminary STAR 
data indicates that b increases when 6 grows, for a fixed Jr, which is in agreement with our qual- 
itative expectations based on the discussion in sec. H] Also consistency between the STAR [1] and 



PHENIX 4l|, |43] measurements is yet not established, as they have both a different experimental 



coverage and analysis method. 
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A. MC Glauber implementation 

We do not use an established Monte Carlo Glauber code, but a very simple implementation 
of our own. For each configuration we first generate two configurations of nucleons according to 
a Woods-Saxon distribution with i?^ = 6.38 fm and surface diffuseness d = 0.535 fm. We then 
draw randomly an impact parameter vector b±. If a nucleon is at a transverse distance of less than 
\/ (^NN /t^ with umn = 41 mb of a nucleon in the other nucleus, it is considered as a participant. 

Once we have a value A'^p, one then has to generate random variables according to the distribu- 
tion eq. ()44|) . For this purpose one must diagonalize the correlation matrix ()44p . This is effectively 
done in the following way. One generates 3 independent Gaussian random numbers ^, ^± with zero 
mean and the variances 

({'> = 1 (46) 

& - (47) 



One then solves Ap b r from 



Af-Ab = V2\/N + aN C (48) 
Af + Ab = Jn + {a + 2K)N^ + (49) 



Aij = V^A^R + (a + K)Nn^ (e+ - e-) • (50) 

It is straightforward to verify that this procedure gives variables N-p,Nb,N^^ with the desired 
correlations, eq. ([illl . 
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